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This contribution briefly reviews some recent work demonstrating the partial break- 
down of the colloidal fluid <-*■ atomic fluid analogy. The success of liquid state the- 
ory for atomic fluids stems in part from the van der Waals picture, where steric 
interactions dominate the structure, and attractive interactions can be added as a 
perturbation. For complex fluids described by effective potentials, this picture may 
break down. In the first example discussed, depletion potentials in non-additive 
hard-sphere mixtures are shown to be surprisingly complex, leading to fluid struc- 
ture and fluid-solid transitions dominated by properties of the attractive potentials 
instead of by the hard-cores. Many colloidal suspensions, and possibly globular 
proteins, fall into this energetic fluid category. In the second example, the coarse- 
graining of polymers leads to soft-core effective potentials and associated mean field 
fluid behaviour distinguished by a breakdown of the virial expansion, an equation 
of state that is nevertheless nearly linear in density, and correlation functions well 
described by the random phase approximation. 
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1. Introduction 

Integrating out a subset of the degrees of freedom (i.e. coarse-graining) is the flrst 
step in many analyses of soft matter systems. For colloidal and polymeric sus- 
pensions, this procedure often leads to effective potentials, and much progress has 
emerged from exploiting the analogy between these potentials and the potentials of 
atomic and molecular fluids. In fact, the basic philosophy behind effective potentials 
is that the initial effort in deriving them is recouped when they are input into the 
well-oiled machinery of liquid state theory. However, coarse-graining on the wide 
range of length scales available in soft-matter systems leads to a much richer class 
of potentials than those found for their atomic and molecular counterparts (Likos 
2001). Since liquid state theory was originally derived and optimized for the latter, 
this immediately implies the possible breakdown of the soft-matter atomic fluids 
analogy. 

This paper will focus on a careful derivation of the effective pair potentials 
for colloid-colloid, polymer-colloid and pure polymer suspensions, with a special 
emphasis on examples where intuition gleaned from the atomic fluid analogy begins 
to fail. 
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In the first example discussed, depletion potentials are derived for highly asym- 
metric non-additive binary hard sphere (HS) mixtures. It is argued that these po- 
tentials can describe a much wider class of asymmetric binary mixtures and lead, in 
the low density regime of the larger species, to energetic fluids, where the structure 
and crystallization behaviour is dominated by the effective potentials rather than 
the hard cores. For this reason, a careful derivation of the effective potentials is 
particularly important. 

The second example concerns modelling polymers as "soft colloids" . These are 
represented by potentials without hard cores, leading to the concept of mean-field 
fluids, with behaviour quite different from the HS paradigm that underpins the 
theory of simple liquids. 

2. The van der Waals picture of fluids 

But first, let us take a careful look at the theory of simple liquids (Hansen & 
McDonald 1986): Why has it been so successful? Perhaps the primary factor is the 
surprisingly widespread applicability of the HS fluid as a model of steric effects in 
atomic and molecular fluids. Attractive forces can then be added as a perturbation 
to a HS reference system. This approach is sometimes called the van der Waals 
picture of fluidsf , and a clear summary statement appeared in an influential review 
(Chandler et al. 1983): 

According to the van der Waals picture, the average relative arrange- 
ments and motions of molecules in a liquid (that is the intermolecu- 
lar structure and correlations) are determined primarily by the local 
packing and steric effects produced by the short-ranged repulsive inter- 
molecular forces. Attractive forces, dipole-dipole interactions, and other 
slowly varying interactions, all play a minor role in the structure, and 
in the simplest approximation their effect can be treated in terms of a 
mean-field - a spatially uniform background potential - which exerts 
no intermolecular force and hence has no effect on the structure or dy- 
namics, but merely provides the cohesive energy that makes the system 
stable at a particular density or pressure. 

Historically, the 1957 discovery of a first order freezing transition in a pure hard- 
sphere system (Alder & Wainwright 1957; Wood & Jacobson 1957) really got the 
metaphorical ball rolling, but other early highlights include: 

• Freezing: Longuet-Higgins & Widom (1964) and Widom (1967) showed that 
attractive interactions only mildly perturb the HS freezing transition near the 
triple point. 

• Structure: Ashcroft and Lekner (1966), and Verlet (1967), modelled the 
structure factor S{k) of liquid metals and many other fiuids by the HS S{k) 
at an appropriate effective HS diameter. Hansen and Verlet (1969) then con- 
nected structure to freezing by deriving a criterion — the first maximum peak 
of S{k) has a value of about 2.8 near the freezing transition — that holds not 
only for HS fluids, but also for a much wider set of atomic fluids. 

t The term van der Waals plus the term fluid appear in various combinations in the literature. 
Its use here follows the dcflnition found in the Chandler, Weeks & Andersen (1983) review. 
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These ideas were put on a firmer footing by the Barker-Henderson perturbation 
theory (Barker & Henderson 1967), and the more systematic Weeks-Chandler- 
Andersen theory of hquids (Chandler & Weeks 1970; Weeks et al. 1971). Both 
theories provide ways of performing quantitative calculations for fluids based on a 
HS reference system plus a perturbative term arising from the attractions. That 
these concepts can also be extended to molecular fluids with non-spherical hard- 
core reference systems is illustrated by a second, more whimsical, citation from the 
Chandler, Weeks, & Andersen (1983) review: 

Similarly the arrangements of molecules in liquid benzene are similar to 
the average arrangements of neighboring Cheerios in a bowl of breakfast 
cereal, and a solution of argon in benzene should be similar to the 
structure achieved when blueberries are mixed with Cheerios. 

The upshot of all this is that near the triple point, the region for which liq- 
uid state theory was optimised, the structure and thermodynamics of atomic and 
molecular fluids is dominated by the underlying hard-core system; the effect of 
attractive potentials is quantitative^ but not qualitative (at least for freezing and 
structure). This is no longer true for complex fluids. Instead, as we shall see in the 
following sections of this paper, effective potentials may induce qualitative changes 
in the structure and thermodynamics of soft matter systems. The van der Waals 
picture breaks down. 

3. Example A: Depletion potentials lead to "Energetic" 

Fluids 

(a) Depletion potentials for non-additive binary hard-sphere mixtures 

Depletion interactions, the effective entropic potentials induced between the 
remaining particles when one (repulsive and typically smaller) component is inte- 
grated out, were first described by Asakura & Oosawa (1954, 1958). They remained 
largely unexplored until around 25 years ago when the experiments of Vincent and 
co-workers (Li-In-On et al. 1975) and the theoretical work of Vrij (1976) rekindled 
interest. Since then, there has been a steadily increasing body of work exploiting 
the "tunability" of depletion potentials: the well-depth is typically proportional to 
the osmotic pressure of the small particles while the range is related to their di- 
ameters, leading to rich and interesting experimental phenomena (see e.g. Crocker 
et al. 1999; Poon et al. 1999; Rudhart et al. 1998; Verma et al. 1998) To study 
this further, consider a binary HS mixture of large (species 1) and small (species 2) 
spheres. The complete description of a binary hard-sphere model demands not only 
the specification of two hard-sphere diameters, cti and (72, but also the specification 
of the cross-diameter (T12 denoting the distance of closest approach between two 
dissimilar spheres. This is traditionally written as: 

ai2 - i(ai+a2) (1 + A) (3.1) 

When A = 0, the cross-diameter is simply the sum of the two radii, exactly what 
one would expect on purely geometric grounds; the model is termed additive, and 
follows the traditional Lorcntz mixing rule (Hansen & McDonald 1986). In nature. 
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however, systems will rarely be exactly Lorentz additive, and in fact positive (A > 
0) or negative (A < 0) non-additivity will be the rule rather than the exception. 
This is further illustrated in figure 1. Each large particle excludes a volume v = 
TTai2/Q from the small particles. When two large particles approach to a distance 
less than 2(Ti2, some of this volume is doubly excluded, and the small particles 
can gain free-volume, and therefore entropy. Integrating out the small particles 
translates this entropy gain into an interaction between the large particles, the 
so-called depletion potential. If 2h > cr2, i.e. the effective exclusion diameter of 
small particles near a large particle is larger than their mutual exclusion diameters, 
then A > and the system displays positive non-additivity. Conversely, if 2h < 
02 then A < and the system displays negative non-additivity. Experimental 
systems which display negative non-additivity include sterically or electrostatically 
stabilized binary colloid mixtures (Louis et al. 2000a). Polymer-colloid mixtures 
typically show positive non-additivity, and in the extreme limit, (T2 — > with h 
finite, the depletion potential {w^'^\Rij) of the appendix) reduces to the Asakura- 
Ooswawa form (Asakura & Oosawa 1958, Vrij 1976): 

^^-('^^ = -^^f (--)^ - + 1^ {£)] ^'-'^ 

in the range ui < r < 2(7i2; here p2 = N2/V. The value at contact is given by: 

PVAoir = ai) = -P2J (^cT,{2hf + ^i^hf^ , (3.3) 

which docs not depend explicitly on the small particle diameter 02 or the packing 
fraction 772- 

For finite 02 the effective pair potentials reduce to (3.2) and (3.3) only in the 
/92 ^ limit. Deriving quantitatively accurate depletion potentials at finite 02 
and p2 for non-additive systems has become possible due to some important new 
developments (see e.g. Roth et al. 2000a). Briefly, the method works like this: 
for a given small particle fugacity 22, the pair-contribution to the exact effective 
potential, u;(^^(r; Z2), is given by the difference in grand potential between a system 
with two large-spheres at a distance r, and the same system with the two large 
spheres at r = 00. Within a density functional theory (DFT) approach, this can be 
rewritten as: 

Pw^^\r;z2)^ lim (c[^\oo) ~ c[^\r)) , (3.4) 

where c"i\r) = —pSF^xlpi, p2]/Spi, Fex[pi, P2] being the excess intrinsic free-energy 
functional of the binary mixture. By employing the Rosenfeld fundamental measure 
theory DFT (Rosenfeld 1989), quantitatively accurate depletion potentials can be 
derived for both additive (Roth et al. 2000a) and non-additive mixtures (Roth & 
Evans 20016). 

The Roth DFT method extracts the effective pair-potential from the one-body 
correlations as shown in equation (3.4). On the other hand, direct functional differ- 
entiation of the same Rosenfeld DFT leads to two-body correlations that are equiva- 
lent to the Percus-Yevick (PY) approximation in the homogeneous limit (Rosenfeld 
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1989). But, as shown in the appendix, PY results in a poor representation of the 
effective pair potentials. Curiously, the same DFT approach provides quantitatively 
accurate effective potentials by one route, and rather poor potentials by another 
route. For this reason, some care must be taken when choosing a particular route 
to thermodynamics or effective potentials from a given (approximate) DFT. 

The effect of non-additivity on the depletion potentials is shown in figure 2(a), 
where the parameters are chosen such that p2h^ is kept constant, but U2 is varied. 
Note in particular that for positive non-additivity (tT2 < 2h) the contact value 
changes very little. In contrast, for negative non-additivity {a^ < 2/i), the contact 
value increases markedly and can even become positive. The non-additivity can also 
be varied by keeping the small-particle diameter a2 and the packing fraction 772 = 
^^92(^2!^ fixed, as shown in figure 2(b). For positive and negative non-additivity, the 
oscillations due to the solvation shells remain more or less the same, but the well 
depth changes markedly. These trends can be understood even at the very simple 
level of equation (3.3), valid only as p — *■ 0, since changing A while keeping (T2 fixed 
corresponds to changing h and therefore the contact value. 

Non-additivity has an even more dramatic effect on the virial coefficients. The 
additive case seems to be marginal since a very small negative or positive non- 
additivity markedly changes the behaviour of B2 away from the additive value 
(Louis & Roth 2000e). A connection to phase behaviour can be made through 
the recent observation of Vliegenthart and Lekkerkerker (2000), who showed that 
32/82^ w — 1.5 near the critical point of a wide variety of fluid systems. The large 
effect of non-additivity on the virial coefficients found by Louis and Roth (2000e) 
can therefore rationalize the large effect of non-additivity found in previous direct 
studies of the fluid-fluid spinodal line (Biben & Hansen 1997, Dijkstra 1998, Louis 
et al. 2000a). Positive non-additivity strongly favours phase-separation, while even 
a very small negative non-additivity has the opposite effect. 

More realistic two-component systems include attractive or repulsive potentials 
Vij (r) in addition to the hard-core steric repulsion. These can be mapped onto the 
non-additivity in the following way: An attractive cross-term Vi2{r) or a repulsive 
small-small interaction V22{i') corresponds to A > 0, while a repulsive Vi2{r) or 
attractive V22(i') lead to A < 0. In this way the large-small or small-small inter- 
actions can be used to "engineer" a very wide variety of effective potential shapes 
and associated fluid behaviours (Louis & Roth 2000e). 

{h) Energetic fluids: Structure and thermodynamics 

(i) Thermodynamics 

Effective depletion potentials can have a much shorter range than the hard-core 
diameter of the larger colloidal species. This leads to perhaps the best-known break- 
down of the simple atomic fluid analogy, namely the metastability of the fluid-fluid 
transition w.r.t. the fluid-solid transition for short-range potentials (Cast et al. 
1983; Hagen & Frenkel 1994). As demonstrated in flgure 3, this leads to crystalliza- 
tion at much smaller large-sphere packing fractions than found for the archetypical 
freezing transition in one-component HS fluids. The relative insensitivity of the 
crystallization line to the form and range of the potentials shown in flgure 3 is 
particularly striking. The determining factor is mainly the depth of the attractive 
well at r = rmin, from which, for the range of potentials probed here, an approxi- 

Article submitted to Royal Society 



6 



A. A. Louis 



mate crystallization criterion can be derived (Louis et al. 2000a): the liquidus line 
broadens to about 50% of the packing fraction at freezing for pure HS's when 
/3K//(w)«2.4±0.3t 

The insensitivity of crystallization to the potential range can be understood from 
the trends in the free-energy curve. With increasing well-depth, the fluid branch of 
the free-energy in figure 4 is only mildly affected, while the solid branch develops a 
deep minimum, driven by close contact of the potential wells at a very high packing 
fractionf . Equilibrium between a very dense solid and a dilute fluid can now easily be 
achieved. (This can be seen in figure 4 by using the common tangent construction.) 
|. Therefore, when the liquidus line broadens, crystallization is driven primarily 
by changes in the close-packed solid branch of the free energy, which probes well 
depths at the potential minima, but is largely independent of other features of the 
potentials. This is the origin of the semi-universal crystallization law found by Louis 
et al. (2000a). 

(ii) Structure 

In the regime where the crystallization curve broadens to lower and lower pack- 
ing fractions of the large-particles, the structure can be well approximated by the 
very simple form: g{r) = exp[— /Jw'^' (r; Z2)] (Louis 2000fo), where w^^^ (r; 2:2) is the 
two-body contribution of the exact effective potential (EEP) described in the ap- 
pendix. As demonstrated in figure 5, this works remarkably well, even for densities 
as high as rji = 0.25. Because the fluid-fluid transition is metastable and crystal- 
lization typically sets in only for a well-depth of order (3V{rmin) ~ 2.4, the g{r) 
near contact can easily approach values greater than 10 in the fluid phase, which is 
quite different from the value for a HS reference system at the same overall density. 

If experiments could directly access g{r) , this would lead to a very simple method 
to extract the effective potentials. However, experiments on pair-correlations typi- 
cally access S{k), where the effect of the potentials is much less clear cut, as demon- 
strated in figure 6. In contrast to the van der Waals picture, the effective potentials 
alter the S{k) quite markedly from the HS reference system, but similarly to the van 
der Waals picture, where deriving the attractive contribution to the potential from 
an inversion of S{k) is very difficult (Reatto 1986), quite a number of different effec- 
tive potentials may result in similar behaviour for S{k). The potentials in figure 6 
were chosen to have nearly equal virial coefficients B2- Since 2B2 = — /(O), where 
f{k) is the Fourier transform (FT) of the Mayer function /(r) = exp[/3ii(r)] — 1, 
(which in turn provides a good approximation for the total correlation function 

t As the potential becomes more short-ranged, the trend is towards somewhat higher values of 
PV^f f{rmin)- Nevertheless, this very simple criterion describes the dominant effect in the range 
of most interest to experiment. 

t This well in the crystal free-energy branch also drives the solid-solid transition 
t While the solid must be treated accurately, small errors in the fluid-free energy have little 
effect on the fluid-solid transition in the energetic fluid regime, in marked contrast to the pure 
HS case where both branches of the free energy must be treated accurately. The free-energies 
in figure 4 were generated with first order thermodynamic perturbation theory. They reproduce 
the fiuid-solid phase-behaviour of short-range potentials rather well, but fail miserably for the 
(metastable) fiuid-fluid transition (Dijkstra et al. f9996, Louis 2000a). First order perturbation 
theory even fails to correctly describe the second virial coefficient, so it is not surprizing that it 
breaks down for the fluid-fluid line. 
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h(r) = g{r) — 1 in this regime) it is not surprizing that S{k) = 1 + ph{k) should be 
similar for different potentials as long as B2 is held constant.^. 

The fluid-solid transition is mainly determined by the potential minimum, and 
is therefore not directly related to B2 (compare, for example, the virial-coefRcients 
of the potentials depicted in the inset of figure 3). This implies that S{k), which is 
mainly determined by will vary significantly for different potentials along the 
fluid-solid line. Therefore the Hansen- Verlet criterion, or any other similar criterion 
based on S{k), will not hold in this "energetic fluid" regime. 

(iii) The energetic fluid picture vs. the van der Waals picture of fluids 

The term energetic fluid applies to the fluid phase at low overall packing fraction 
of the large-particles (say rji < 0.25), as shown schematically in figure 3. This regime 
opens up in short-range potential fluids when the fluid-fluid phase-line becomes 
mctastablc to the fluid-solid line, which then broadens out to very low packing 
fractions. While it is hard to flnd atomic or molecular fluids with a potential range 
short enough to qualify for this nomenclature, many solutes in solution arc governed 
by relatively short-range attractive potentials. Through McMillan-Mayer solution 
theory (McMillan and Mayer f945) they can be viewed as effective atomic fluids. 
But in contrast to the van der Waals picture described in section 2, the attractive 
potentials of such energetic fluids qualitatively affect both freezing and structure: 

• Freezing: The liquidus line is determined primarily by the well-depth of the 
potential, and less so by other details such as its range or shape. While in the 
van der Waals picture freezing is driven by cntropic excluded volume effects, 
here it is driven by (energetic) properties of the potentials. 

• Structure: The real-space structure of energetic fluids can be qualitatively 
described by the simple relation g{r) = exp[— /3z;(r)], while the k-space struc- 
ture S{k) is well described by the Baxter model at the same i?2- In contrast, 
the S{k) of van der Waals fluids typically resembles that of a HS fluid. 

is mainly determined by the value of the second virial coefficient rather than 
by other details of the potentials. 

Examples of energetic fluids include many colloidal suspensions, where direct 
interactions can be of short range, or else where other (smaller) colloids, micelles, 
or polymers can act as depletants, inducing potentials such as those depicted in 
figure 2. 

Globular proteins may also fall into the same class, but their interactions are 
no longer spherically symmetric (Durbin & Feher f 996; Lomakin et al. 1998; Neal 
et al. 1998; Piazza 2000; Sear 1999). The radial distribution function and structure 
factor should still display behaviour similar to other energetic fluids, but now with 
an implicit orientational average over the anisotropic potentials. 

Proteins typically crystallize at very low packing fractions, and some show evi- 
dence of a phase diagram topology similar to the one depicted in figure 3 (Frenkel 
2000; Piazza 2000; Rosenbaum 1996; Sear 1999), with a metastable fluid-fluid line. 
There is also evidence that potential contacts in the crystal are very important 

^ This simple picture breaks down near the spinodal line where S{k = 0) is enhanced by 
collective density fluctuations. 
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for protein crystallization (Durbin & Feher 1996), which is consistent with ener- 
getic fluid behaviour and in sharp contrast to what would be expected from the 
van der Waals picture. The semi-universality of crystallization found for spherically 
symmetric potentials depends on similar behaviour of the solid branch of the free- 
energy for different potential shapes. Unfortunately, this can no longer be expected 
for proteins, since (a) the proteins are no longer spherical and (b) the potentials 
are anisotropic (or "patchy"), leading to many possible crystal structures. In fact, 
as emphasized by Sear (1999), conditions (a) and (b) suggest that some proteins 
may be very difficult to crystallize precisely because their potential contacts are not 
commensurate with an allowed crystal structure. So, whether globular proteins can 
be usefully classified as energetic fluids remains to be seen, but they can certainly 
not be classified within the van der Waals picture of fluids. 

4. Example B: Polymers as soft colloids: Mean field fluids 



The approach outlined here, described in more detail in (Louis et al. 2000c, 
Bolhuis et al. 2001), is to coarse-grain the polymers by integrating out monomeric 
degrees of freedom, resulting in a description based on the polymer centres of mass 
(CM). In principle, for a set of A'^i polymers with L monomers each, one could 
calculate the EEP which depends explicitly on the configuration {Ri} of the centres 
of mass. The only differences with the binary-mixture case described in the appendix 
are: (1) there is no explicit additional dependence on the monomer density (or 
fugacity), since that is fixed by the number of polymers in the set {Ri}, and (2) 
there is no direct interaction V{{Ri}); the polymer-polymer effective potential arises 
solely from the coarse-graining procedure. The EEP now takes the form: 



and again, each term could in principle be calculated separately, but this rapidly 
becomes intractable. Instead, one can follow the second route described in the 
appendix and include all higher order terms in equation (4.1) through the density 
dependence of an effective pair potential z/;(r; pi), which exactly reproduces the CM 
pair-correlations at a given density pi ~ Ni/V. To do so, we first perform Monte 
Carlo (MC) simulations of self avoiding walk L = 500 chains on a simple cubic 
lattice and extract the radial-distribution functions g{r) between the CM. These 
are then inverted through an Ornstein-Zernike procedure using the hypernetted- 
chain closure (HNC), which has been shown to be quasi-exact for the potentials 
under consideration (Louis et al. 20006, 2000(i, Bolhuis et al. 2001). Figure 7 shows 
both the radial-distribution functions and the effective potentials inverted from the 
g{r). According to the theorems discussed in the appendix, the effective-potentials 
should, through the pair-correlations and equation (A 9), reproduce the contribu- 
tions to the (osmotic) compressibility induced by the EEP (4.1). Figure 8 shows 
that the total compressibility is very well approximated in this approach, implying 
that, in contrast to the binary fluid case, the volume term contribution (defined in 
the appendix) is small. The difference in volume terms is most likely due to the 
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fact that the number of monomers is fixed by the number of polymers, instead of 
being free to vary as is the case for the small particles in the binary mixture. In 
fact, if the polymers were rigid, their CM correlations would exactly determine the 
compressibility (see e.g. Hansen & McDonald 1986 ch 12). 



The effective potentials shown in figure 7 are radically different from the usual 
hard-core plus attraction form found for atomic systems, immediately suggesting 
different fluid behaviour. Already the radial distribution functions in figure 7 appear 
to be quite different from their atomic- fluid counterparts. The structure of fluids 
described by these potentials is well approximated by the simple random phase 
approximation over a surprizingly large density range (Louis et al. 2000(i, Likos et 
al. 2000). Marked differences with atomic fluid behaviour also arise at the level of 
thermodynamics, where the equation of state (EOS) is very well described by the 
following mean-field form; 



As shown in figure 9, this also hold for integrablc potentials which diverge at the 
origin. Equation (4.2) demonstrates that what matters is the value of the FT of 
the potential at the origin, /3u(fc ~ 0), which is finite for all 3 potentials, and, 
in figure 9, chosen to have the same value to ensure that the three EOS are very 
close. The quasi-linear behaviour of the EOS resembles that of an atomic fluid in the 
second virial approximation, but, in fact, the three potentials have virial coefficients 
differing by about 10%. More seriously, the virial expansion has a surprisingly small 
radius of convergence so that adding higher order virial coefficients results in a much 
poorer approximation (Louis et al. 2000d) , as demonstrated by the third order virial 
expansions in figure 9. 

Instead of the CM, end-points or mid-points could also be used to construct an 
effective particle picture of interacting polymer solutions. For example, the mid- 
point representation would be very similar to the two-arm limit of a star-polymer, 
for which a number of results have been recently derived (Likos et al. 1998; Likos 
2001; Watzlawck et al. 1999). In figure 10 the / = 2 limit of the star-polymer 
potential is compared to a more recent expression for the mid-point-mid-point 
interaction (DzubicUa et al. 2000). At first sight the two appear very similar, but 
the slope of their respective EOS differ by a factor of four! The reason for this can 
be seen in figure 11: Since (ivik = 0) is proportional to the integral over r'^(3v{r), 
it is the small differences (much less than fc^T) in the tails of the two (3v{r) that 
determine the large differences in the EOS. Interestingly the two radial-distribution 
functions are again quite similar (although the associated structure factors are not). 



A necessary, but not sufficient condition for mean field fluid (MFF) behaviour 
is that I3v{k = 0) is finite. If /3u(fc = 0) is relatively large, then MFF behaviour, 
with its quasi-linear EOS and correlation functions well-described by the random 
phase approximation, only sets in for large enough densities. At lower densities the 
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structure can still resemble that of a HS fluid, but the topology of the phase-diagram 
is quite different (Lang et al. 2000, Louis et al. 2000(i). Depending on whether the 
FT of the potential oscillates or not, such fluids may show either re-entrant melting 
or else a clustering transition for large enough values of (3v{k = 0) (Likos et al. 
2000). 

Systems with values of f3v{k = 0) significantly below the value that can induce 
freezing or clustering transitions are MFF at all densities. The associated simple 
equation of state and fluid-structure are markedly different from the behaviour of 
fluids in the van der Waals picturef. 

5. Conclusions 

While the van der Waals picture of fluids lies at the basis of many successful liquid 
state theories for atomic and molecular fluids, it breaks down for the richer class of 
potentials arising from coarse-graining procedures in complex fluids. Nevertheless, 
a judicious choice of techniques drawn from the theory of simple liquids may still 
provide insight when combined with a careful derivation of effective potentials. 
Some recent examples of this programme include the phase-diagram of additive 
(Dijkstra et al. 1999a) and non-additive (Louis et al. 2000a) asymmetric binary HS 
mixtures, the structure and phase-behaviour of star-polymers (Likos et. al. 1998; 
Watzlawek et. al. 1999), and star-polymer colloid mixtures (Dzubiella et. al. 2000), 
and the structure and phase-behaviour of a pure-polymer system (Louis et. al. 
2000c; Bolhuis et. al. 2001). In each case the effective potential picture was the key 
to clarifying the underlying physics. 

Similarly, in this contribution a derivation of the effective potentials in binary 
hard-sphere fluids demonstrated the key role of non-additivity in determining the 
shape of the effective depletion potentials, and the associated phase behaviour. 
When such potentials are short ranged, they open up a region in the phase-diagram 
where the fluid behaviour deviates significantly from the van der Waals picture, so 
much so in fact, that a new nomenclature, energetic fluids, has been introduced. In 
the energetic fluid regime, the structure can be described by a very simple form: in 
real space g[r) = exp[— /Jw*-^^ (r; p2)\, while in reciprocal space S{k) is largely deter- 
mined by the virial coefficient B2. For a wide variety of potential shapes, the liquidus 
line broadens to low large-particle densities when well-depth (iw^'^\rmin] P2) ~ 2.4. 

When linear or star polymers in solution are represented as "soft colloids" cen- 
tred around their mid-points or their CM, the resultant picture leads to effective 
potentials with a finite value of /3v{k = 0). Again, such fluids do not follow the van 
der Waals picture. For example, their EOS closely follows a mean-field linear form 
(3P/p w 1 4- \/2pPv{k = 0), and their structure is well described by the simple 
random phase approximation closure: behaviour best classified under the moniker 
mean field fluids. 

Most of this work was done in close collaboration with J. P. Hansen, P.G. Bolhuis, R. 
Roth, R. Finken, and E.J. Meijer, and has appeared in our joint publications listed in the 
references. I thank David Rowan for a critical reading of the manuscript and gratefully 
acknowledge financial support from the Isaac Newton Trust, Cambridge. 

t Perhaps I'm beating a straw man here, but the distinctions are hopefully helpful nonetheless 
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Appendix A. Two ways of deriving effective potentials for 

binary mixtures 

Consider a binary mixture of large (species 1) and small (species 2) spherical par- 
ticles. Integrating out the smaller component to derive a new one-component fluid 
interacting through effective depletion interactions is a useful way to treat such 
systems. This appendix will focus on two popular ways to derive these effective 
potentials: 

(i) Method 1: Exact effective potential 

For binary mixtures the exact effective potential (EEP) is most easily analyzed 
in the semi-grand ensemble, where the number of large particles A^i and the fugacity 
Z2 of the small particles is fixed (Lekkerkerker et al. 1992, McMillan & Mayer 
1945)t. Given a set of Ni large particles fixed at positions {Ri} in a volume V, 
the small particles are integrated out by calculating their partition function in the 
fixed external field generated by the large particles. This results in an effective 
grand potential for the small particles of the form: 

n{N^,z2,V-{R,}), (Al) 

which depends parametrically on the large-particle positions {RfW- To make further 
progress we follow the analysis of Dijkstra et. al (1999a, 1999&, 2000), and decom- 
pose this grand potential into n-body terms. Using this mapping, a one-component 
system can be derived with an an effective interaction of the form: 

vii{{Ri}) + n{Ni,z2,V- {Ri}) - ^("'(^^1, ^2, V) + W'ff{{Ri}- Z2), (A 2) 

where vii({i?i}) is the direct interaction between the large spheres that is already 
present in the original two component system, y(°)(A^i, Z2, V) is the so-called vol- 
ume term (which depends on the set {-Ri} only through the total number of large 
particles iVi, but is independent of their relative positions)^, and W^^ ^ {{Ri}] Z2) 
is the EEP which can be further expanded as: 

W^ff{{R};z2)=Y. {''i^{R^j)+w(^\R,f,Z2)) + ^ (v''^ {R^^u) ^ w^^^\R,,k; z^)) + . . . 

i<j i<j<k 

(A3) 

The pair term is simply the sum of the pair contribution to the direct interaction, 
v'^i{Rij), and the effective interaction w^'^\Rij] Z2) induced by integrating out the 
small spheres. It is precisely this pair term which is measured in the laser tweezer 
experiments (Verma et al. 1998, Crocker et al. 1999). The triplet term is a similar 
sum of direct and induced terms with Rijk the standard 3-body coordinates and 

f Charged systems arc best analyzed in the the canonical ensemble. An early example would 
be the effective potentials in liquid metals (Ashcroft & Stroud 1978) 
X Explicit temperature dependence is omitted in this section 

^ For binary fluids the volume term is given by: V^°\Ni,Z2,V) = -P2(Z2)V + NiwW(z2) 
where p2 is the pressure of a pure component 2 system at fugacity Z2, and w^^\z2) is the grand 
potential difference due to adding one large component 1 particle to the same system. One body 
terms which depend on the configuration {-Ri} are zero duo to translational invarianco. 
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so forth for n-body interactions. (Since the individual terms w^^^Rij,,,) depend 
only on the large-small and small-small interactions, their form is independent of 
which need not be decomposable as a sum of pair terms.) The induced 
pair term can be rigorously defined as the difference between the grand potential 
of equation (A 1) for two particles a distance Rij = |Ri — Rj| apart, and the grand 
potential when the two particles are an infinite distance apart. From (A3) this 
reduces to w^^\R^j]Z2) = n[Ni = 2,Z2,V;R,j) - V^°^{Ni = 2,Z2,V). Similarly 
the induced three-body potential can be defined as: 

w^^\R,jk;z2) = niNi=3,Z2,V;R,jk)~V'-"HNi=3,Z2,V) 

- w'-^\R,y,Z2)-w^^HR,k;z2)-w'^^\Rjk;z2), (A 4) 

i.e. it is that part of the interaction induced by 3 large spheres which cannot be 
described by volume and pair interaction terms alone. By continuing in similar 
fashion for higher and higher order terms the full EEP can be built up such that 
the exact free energy of the mixture is given by: 

F{Ni,z2,V) = V^°\m,z2,V)+F-ff{Ni,Z2,V), (A5) 

where F'^^-f {Ni, Z2,V) is defined as the free energy of a one-component system 
interacting through the EEP, i.e. exp[~(3F'^^-f] = Tri cxp[— /^VF'^-'''^]. In other words, 
by exactly integrating out the small particles for an arbitrary configuration {Ri\, 
the original two-component partition sum has been rewritten as a weighted sum 
over large particle configurations only. The EEP describes the weighting of each 
individual large-particle configuration in the effective one-component partition sum, 
while the volume term adds the left over contributions independent of the large- 
particle configurations. 

Calculating the EEP to all orders is usually impractical. Instead what is often 
done is to truncate the series (A 3) and retain only the pair potential. In many 
cases this is not such a bad approximation, and for a few systems it is even exact. 
For example, the EEP for the Asakura Oosawa (AO) model (Asakura & Oosawa 
1958, Vrij 1976) with size-ratio q = 02! o\ < 2/\/3 — 1 — 0.1547 is exactly described 
by the pair termf. Similarly, for 0.1547 < g < ^/3/2 - 1 ~ 0.2247 the EEP is 
exactly described by the pair term and a triplet term. But even for this rather 
simple model, the exact form of u>^^^ (Rijk', Z2) is very tedious to calculate (Goulding 
2000). Thankfully, simulations have shown that neglecting the higher order terms 
for the AO model often works admirably well even for size ratios as large as g = 0.5, 
(Meijer & Frenkel 1994; Dijkstra et al. 1999a, 19996). 

(ii) Method 2: Inverting the pair correlations 

The correlations functions in an effective one-component system interacting 
through the full EEP (A 3) arc equal to the correlation functions between the large 
particles in the original two-component system, as long as both are at the same 
state point (see e.g. Dijkstra et al. 2000)|. For example, S{k) = Sii{k), where the 

t Another oven simpler example with only a pair term is the exactly solvable lattice model of 
Frenkel and Louis (1992). 

X The volume terms do not directly contribute to the pair and higher order correlations al- 
though they may contribute indirectly through changing the phase-behaviour of the system. For 
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first structure factor is that of the effective one-component system and the second 
is that of the original two-component system. This equivalence also implies that 

lim gn(r;z2) =cxp [-/3f«fi^(r) +«;(2)(r;z2))l ' (^6) 
Pi^O L \ /J 

where gii(r]Z2) is the radial distribution function of the large-spheres, f. Similar 
expressions can be derived to link higher order correlation functions to higher order 
terms in the EEP (A 3). Hence, there is a direct link between the EEP of the effective 
one-component system and the pi ^ limit of the fluid correlation functions in 
the homogeneous phase of the original two-component model| 

The relation (A 6) may be generalized to finite densities of species 1, because 
there exists a one-to-one mapping: 

g{r] p) ^ w{r] p) (A 7) 

between a given pair distribution function g{r) at density p and a unique two-body 
pair potential w{r] p) which reproduces g{r) irrespective of the underlying many- 
body interactions in the system (Henderson 1974; Chayes et al. 1984). Since at 
finite densities the radial distribution function g(r;pi,Z2) = guir; pi, Z2) includes 
contributions not only from the pair-interactions, but also from higher order con- 
tributions to W^f^{r; Z2, {Ri}), the effective pair-potential includes these terms in 
an averaged way. It can be written as: 

w(r;pi,Z2) = wii(r) +w{r;pi,Z2) (A 8) 

which is connected to equation (A6) through lim p-^ ^0 w{r; pi, Z2) = w^^'(r;z2). 

The price paid for including the effect of all higher order terms is to introduce a 
density dependence in the pair-potential, but the payoff is that the pair-correlations 
are exactly reproduced (but not the triplet or higher order correlations). Thermo- 
dynamics can then be extracted through the compressibility relation: 

V 9pi J ''^^ S{k) 

Note, however, that volume terms can also contribute to the total EOS (Dijkstra et 
al. 2000) , so that equation (A 9) describes the osmotic compressibility due to species 
1 only. In some cases this is only a small fraction of the total compressibility of the 
full two-component system (Louis et al. 1999; Dijkstra et al. 2000). 

How does one perform the g{r] p) <-!■ w{r; p) inversion? If the full EEP includes 
only pair terms (as is the case for the AO model with q < 0.1547), then w{r; pi, Z2) = 

the binary mixtures considered here they don't affect phase-boundaries (Dijkstra et al. 1999a), 
but for other systems they can (van Roij & Hansen 1997; van Roij et al. 1999; Graf & Lowen 
1998; Warren 2000), and so must be taken into account to correctly describe the equivalent state 
points of the original two-component and effective one-component systems 
t Most information about gi2{k) and 922 C^) is lost in the mapping. 

I Note that equation A 6 is different from the potential of mean force which is typically defined 
as f3wP"^f (r) = — ln[g(r)] for any density and is strictly speaking not a pair potential, but a 
restatement of the pair-distribution function. For example, a system directly interacting through 
a l3wP"^f (r) derived at finite p will not have the same correlations as the original system from 
which the potential of mean force was derived. 
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w^^^r, Z2) at all densities. If it only includes pair and triplet terms (as for the AO 
model with 0.1547 < q < 0.2247) then to a good approximation (Reatto & Tau 
1987; Attard 1992): 

w{ri2:Pi,Z2) - w^^\n2;z2) ^pij drs [c-'"'''''^^''^^''-^--^) - l] g(ri3).g(r23). 

(A 10) 

More generally, one needs (i) a method to generate the exact 511 (r) for the two 
component system and (ii) an inversion method to extract j3v{r) from g{r). Inver- 
sion methods based on the Ornstein-Zernike relations exist (Reatto 1986, Zerah 
and Hansen 1986), but these are very sensitive to the underlying approximations 
and the quality of the original pair-correlations used as input. 

As an example of the difficulties involved in (i), consider the popular Per- 
cus Yevick (PY) approximation which is exactly soluble for binary HS mixtures 
(Lcbowitz & Rowlinson 1965). On the one hand PY approximates the EOS very 
well, but on the other hand the PY approximation to the large-particle correlation 
function in the pi,p2 ^ limit reduces to: 

gP^{r)^\-(3VAo{r) (All) 

Pl,P-2—'0 

instead of the correct exponential form: exp[— /3V}io(?')] (note that 'w^'^\r\ P2) ~ 
Vao('') in this limit). For finite p2, w^'^\r; P2) begins to deviate from the AO form, 
but PY still approximately linearizes the exponential, as illustrated in figure 12. In 
spite of the fact that PY successfully describes the EOS, it clearly fails quite badly 
for the effective potential, especially near contactf, suggesting that an inversion of 
PY at finite densities should also fail. More generally, using 2-component integral 
equations to derive effective pair potentials or phase-behaviour in the "colloidal 
limit" (small y, large 772 and small 771) is fraught with difficulty. For example, Biben 
et al. (1996) showed how two self-consistent closures that work extremely well for 
one-component systems, RY and BPGG, predict quite different locations of the 
fluid-fluid spinodal line in binary HS mixtures. A comparison of the effective pair 
potentials ?«^^^(r; P2) calculated for each closure by these authors can, in fact, ratio- 
nalize the difference, since these effective potentials are what primarily determines 
the phase-behaviour. This suggests that for two-component integral equations in the 
colloidal limit, it is much better to compare the performance for the large particle 
osmotic pressure rather than the performance for the full EOS, which is typically 
dominated by the small-particle contribution. 
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Figure 1. The centres of the small spheres can only approach to within a distance 
h = (Ji2 — ai/2 = ^ai{A + y + Ay) from the surface of the large spheres {y = cr2/o"i)- 
For an additive system a2 = h, but here CT2 < so that the system exhibits positive 
non-additivity. 
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Figure 2. The effect of non-additivity on depletion pair potentials at a 
size-ratio y = a^jox — 0.1. In plot (a) the effective packing fraction 
ri2^^ = 4:/'i-Kp2h? — 0.258045 is kept constant, while in plot (b) the small-particle 
packing fraction 772 = 7r/6p20"2 ~ 0.258045 is kept constant. For the AO potential 
at y = 0.1, this packing corresponds to a second virial coefficient 32/82 '^ = —1.5, 
which is near the (metastable) fluid-fluid critical point (Vliegenthart & Lekkerkerker 
2000). The potentials are from Louis et al. (2000a), and ignore oscillations at a range 
r > 2h. Note how the two ways of varying the non-additivity affect the depletion 
potentials differently. 
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jFi gure 3. Typical phase-diagram for short-range depletion potential systems. The 
meta-stable fluid-fluid lines are for depletion potentials with y — 0.2 and different A 
(Louis et al. 2000a). The fluid-solid lines come from the potentials in the inset (each 
with corresponding line styles), here normalized to the same minimum at contact. 
Note the differences in range and shape of the potentials. The area inside the broad 
dashed lines roughly denotes the stable energetic fluid regime; above the fluid-solid 
curves would be the metastable energetic fluid regime. The fluid-solid lines for the 3 
shorter range potential systems was generated with first order perturbation theory, 
while the y = 0.4 AO potential fluid-solid line comes from the simulations of Dijkstra 
et al. (19996) because 1st order perturbation theory incorrectly exhibits a stable 
fluid-fluid phase-transition. 
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Figure 4. Normalized free energies per unit volume for the additive HS poten- 
tial used in (Dijkstra et al. 1998), at a size-ratio y = 0.2. The packings of the 
small-spheres are ri2 = 0.02, 0.08, 0.16, 0.22 and ,0.26, which correspond to well-depths 
of l3V{r = o-i) = 0.183,0.768,1.597,2.224,2.632, respectively. Both branches of the 
free energy curve decrease with increasing potential well-depth, but the effect is much 
more pronounced for the crystal branch. 
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Figure 5. The exponential form g{r) = exp(— /JVaoC?") is a semi-quantitative approx- 
imation for these two state-points in the energetic fluid regime. When rji = 0.25, the 
system is near the fluid-sohd hne. The g{r) were generated with the Percus Yevick 
(PY) approximation, which is quantitatively accurate in this regime (see e.g. Dijkstra 
et al. 19996)). 
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Figure 6. The structure-factors S{k) at rji = 0.13 for the potentials in the inset. Note 
that the shape of S{k) is poorly approximated by the pure hard-sphere form, but well 
represented by the Baxter model S{k). The two shorter range potentials have the 
same B2 as the Baxter model, while the two longer range potentials have a B2 about 
7% more negative, which fits to almost the same S(0). The differences in the S{k) 
due to the shape and range of the potentials are not much larger than the resolution 
of the best experiments on colloids in solution. As in the previous figure, the S{k) 
were determined by the accurate PY approximation. 
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Figure 7. MC simulations are used to generate the g{r) for L — 500 SAW poly- 
mer chains on a cubic (240"^) lattice for densities from p ~ (A'^ = 2 polymers) to 
p/p* = 8.7, {N = 6400) Here p* = 4/37ri?g, and p/p* < 1 denotes the dilute and 
p/p* > 1 the semi-dilute regimes for polymers in a good solvent. These radial dis- 
tribution functions are inverted using an Ornstein-Zernike procedure to obtain the 
effective potentials between the polymer CM. 
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Figure 8. The equation of state derived from the compressibility relation (A 9) accu- 
rately approximates the true EOS measured by direct simulations of the polymers. 
This implies that the volume terms are small. Ignoring the density dependence of 
the potentials by using only the low density form of the pair potential, w{r\ p = 0), 
equivalent to the pair-contribution to the EEP, leads to an underestimate of the 
pressure 
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Figure 9. (a): Three potentials PV{r) which all result in a mean field fluid, (b): The 
EOS Z = PP/ p, here generated by the quantitatively accurate HNC approximation, 
are all very close to the Zmf (equation (4.2)) form (circles). However they are badly 
approximated by a third order virial series Z = 1 + 52/9 + -Bsp'^, here represented 
by squares along each line. In this and the next two plots, line styles in (a) and (b) 
denote corresponding systems. 
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Figure 10. (a) Comparison of two very similar looking potentials that describe the 
midpoint-midpoint interaction, (b) Note how different their EOS are! 




5 10 15 1 2 3 4 



r r 

Figure 11. (a) The differences between the two mid-point potentials in the previous 
figure become much more apparent when they are multiplied by . (b) This does 
not translate into much difference in the g{r) shown here for p — 0.5, but S{k) (not 
depicted here) does change. 
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Figure 12. Comparison of the depletion potential w^^\r; P2) derived from the PY 
approximation and from the quantitatively accurate Roth DFT method discussed in 
the text. PY tends to underestimate the potential strength, leading to a much reduced 
contact value of .gii(r) in the pi ^ limit, as shown in the inset. 
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